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The charged-current quasi-elastic scattering of muon neutrino on a carbon target is 
calculated for neutrino energy up to 2.8 GeV using the relativistic distorted-wave im- 
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■ pulse approximation with relativistic optical potential, which was earlier successfully 



applied to describe electron-nucleus data. We studied both neutrino and electron 
processes and have shown that the reduced exclusive cross section for neutrino and 
^ ' electron scattering are similar. We have also studied nuclear and axial vector mass 

a: _ , 

effects on the shape of Q distribution. The comparison of the (anti) neutrino to- 
tal cross sections per (proton)neutron, calculated for the carbon and oxygen targets 
shows the cross sections for oxygen to be lower than those for carbon. We found 



significant nuclear model dependence of inclusive and total cross sections for energy 
about 1 GeV. 

PACS numbers: 25.30.-c, 25.30.Bf, 25.30.Pt, 13.15.+g 
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c3 ■ I. INTRODUCTION 

The goals of the current and planed set of accelerator-based neutrino experiments [1-7] 
are the precision measurements of the neutrino mass squared difference Am| 3 by measuring 
muon neutrino disappearance, and searching for the last unmeasured leptonic mixing angle 
#13 through the muon to electron neutrino transition. The last oscillation channel, if it 
occurs, opens the possibility of observation matter/anti-matter asymmetries in neutrinos 
and determination of the ordering of the neutrino mass states. The data of these experiments 
will greatly extend the statistics due to extremely intense neutrino beamline. 

To study the neutrino oscillation effects on the terrestrial distance scale, the neutrino 
beams cover the energy range from a few hundred MeV to several GeV. In this energy range, 
the dominant contribution to the neutrino-nucleus cross section comes from the charged- 
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current (CC) quasi-elastic (QE) reactions and resonance production processes. The cross 
section data in this energy range are rather scarce and were taken on the targets, which 
are not used in the neutrino oscillation experiments (i.e. water, iron, lead or plastic). In 
this situation, the statistical uncertainties should be negligible as compared to systematic 
errors in the incident neutrino flux, neutrino interaction model and the detector effects 
on the neutrino events selection and neutrino energy reconstruction. Apparently, these 
uncertainties produce systematic errors in the extraction of oscillation parameters. 

Many experiments try to reduce these uncertainties by using a near detectors. One of 
the option for near detector design is to make the near detector as more segmented and fine- 
grained, using scintillator (carbon) as a target and detecting material. This strategy means 
that one must try to measure the fluxes and cross sections as independently as possible, 
and then use this information to constrain the detector simulation so that the information 
is correctly extrapolated to the far detector. The concern with this strategy is that the 
detector simulation must accurately predicted the detector response. Because the near and 
far detectors are not necessarily of the same target material, a part of the near detector must 
include some of the same target material, so that nuclear effects on the cross sections (QE 
and non-QE) could be taking into account. Among the proposed experiments MINERvA [4] 
and ND280 detector [8] will have the possibility for precise measurements in a wide range 
of energies and for various nuclear targets. 

To model the QE neutrino scattering from a nuclei, the most part of the Monte Carlo 
(MC) event generators [9] are based on the relativistic Fermi gas model (RFGM) [10] with 
Pauli blocking, in which the nucleus is described as a system of quasi-free nucleons with 
a flat nucleon momentum distribution up to the same Fermi momentum pp and nuclear 
binding energy e^. But this model does not take into account the nuclear shall structure, 
the final state interaction (FSI) between the outgoing nucleon and residual nucleus and 
the presence of short-range nucleon-nucleon (AW) correlations, leading to appearance of a 
high-momentum and high-energy component in the nucleon momentum-energy distribution 
in the target. 

The comparison with the high-precision electron scattering data has shown [11] that the 
accuracy of the RFGM prediction becomes poor at low squared four-momentum transfer 
Q 2 , where the nuclear effects are largest. The modern quasi-elastic neutrino scattering data 
(the CC QE event distribution as a function of Q 2 ) [1, 12] also reveal the inadequacies in 
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the present neutrino cross section simulation. The data/MC disagreement shows the data 
deficit in the low-Q 2 (Q 2 < 0.2 (GeV/c) 2 ) region. 

There are many calculations for the QE neutrino charged-current and neutral-current 
scattering cross sections on nucleus, which go beyond the simple RFGM and use more 
realistic description of nuclear dynamics. In calculation of Refs.[13, 14] within the plane- 
wave impulse approximation (PWIA), the short range iViV-correlations were included using 
the description of nuclear dynamics, based on nuclear many body theory. Charged current 
and/or neutral current neutrino-nucleus cross sections were studied within the relativistic 
distorted- wave impulse approximation (RDWIA) in Refs. [15-24], using the relativistic shell 
model approach and taking into account the FSI effects. In Refs. [21, 22] the contribution 
of the short range correlations (SRC) was also considered. The FSI effects were studied 
in Refs. [25-28] within the framework of the random phase approximation, in Refs. [29-31] - 
within a Superscaling approach, and in Ref. [32] - in a GiBUU model. 

In this paper, we calculate the single-nucleon knockout contribution to the exclusive, 
inclusive and total cross sections of charged-current QE (anti) neutrino scattering from 12 C, 
using different approximations (PWIA and RDWIA) and the Fermi gas model. We employ 
the LEA code [33] which was adopted for neutrino reactions. In our approach, the effect of 
the SCR in the carbon ground state is evaluated in the PWIA [34, 35] and the FSI effect 
on the inclusive cross sections in the presence of the iViV-correlations is estimated according 
Ref. [21]. The aims of this work are a) calculation the RDWIA CC QE z/ 12 C cross sections, 
b) investigation of nuclear effects on the Q 2 dependence of the (anti)neutrino cross section, 
and c) comparison of the total cross sections, scaled with the number of neutron/proton in 
the target for (anti) neutrino scattering on the oxygen and carbon targets. 

The outline of this paper is the following. In Sec. II we present briefly the formalism 
for the CC QE scattering process and the RDWIA model. The results are presented and 
discussed in Sec. III. Our conclusions are summarized in Sec. IV. 

II. FORMALISM OF QUASI-ELASTIC SCATTERING AND RDWIA 

We consider electron and neutrino charged-current QE exclusive 



1(h) + A(p A ) - l'(kf) + N( Px ) + B(p B ), 



(1) 
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and inclusive 

l(ki) + A(p A )^l\k f )+X (2) 

scattering off nuclei in one-photon (W-boson) exchange approximation. Here I labels the 
incident lepton [electron or muon (anti) neutrino], and I' represents the scattered lepton 
(electron or muon), hi = (£j, fej) and kf = (ef,kf) are the initial and final lepton momenta, 
Pa = ( £ a,Pa)i an d Pb = ( £ b,Pb) are the initial and final target momenta, p x = (e x ,p x ) is 
the ejectile nucleon momentum, q = (u, q) is the momentum transfer carried by the virtual 
photon (W-boson), and Q 2 = —q 2 = q 2 — u 2 is the photon (W-boson) virtuality. 



A. CC QE neutrino-nucleus cross sections 

In the laboratory frame, the differential cross section for the exclusive electron (cr el ) 
and (anti)neutrino CCQE (<x cc ) scattering, in which only a single discrete state or narrow 
resonance of the target is excited, can be written as 

d 5 a el 



R \Px\ £ x £ f® j (el) w Mel) (n s 

de f dQ f dtt x (2tt) 3 e % Q 4 ^ y ' 

d5 ° cc =R \^M&™ 29 ° L M W *>M (3b ) 



dSfdVtfdVL x (27r) 5 Si 
where Sly is the solid angle for the lepton momentum, Q x is the solid angle for the ejectile 
nucleon momentum, a ~ 1/137 is the fine-structure constant, G ~ 1.16639 xlO" 11 MeV~ 2 
is the Fermi constant, 9c is the Cabbibo angle (cos 6*c ~ 0.9749), L^ u is the lepton tensor 
and wffl and wf& are, respectively, the electromagnetic and weak CC nuclear tensors. 



The recoil factor R is given by 

£ x P x ■ P B 



R — J de x 5(e x + e B - u> - m A ) = 
and e x is the solution to the equation 



£ B P x ■ P x 



(4) 



e x + e B -rn A -uj = 0, (5) 

where e B = ^/ m B + Pbi Pb = Q ~Pxi Px = V £ x ~ m2 i an d m A, t^b, and m are masses of 
the target, recoil nucleus and nucleon, respectively. The missing momentum p m and missing 
energy e m are defined by 

Pm = P x ~ Q (6a) 
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e m = m + m B -m A (6b) 

The leptonic tensor is separated into a symmetrical and an anti-symmetrical components 
that are written as in Ref. [21]. The electromagnetic and weak CC hadronic tensors, 



and wffi are given by bilinear products of the transition matrix elements of the nuclear 
electromagnetic or CC operator jjf 1 ^ ^ between the initial nucleus state \A) and the final 
state \Bf) as 

W^ )(cc) = ^(B f ,p x \jWl"M)(M*™\Bf,Px), (7) 
/ 

where the sum is taken over undetected states. 

In the inclusive reactions (2) only the outgoing lepton is detected, and the differential 
cross sections can be written as 

= ^L^W** (8a) 



defdVtf €iQ 4 

= J- M ^IN^ (8b) 
de f dVtf (2vr) 2 e t 2 ' { } 

where is the inclusive hadronic tensor. The expressions for the exclusive (3) and 

inclusive (8) lepton scattering cross sections in terms of response functions are given in 

Ref. [21]. 

It is also useful to define a reduced cross section 

*- = l^nj^^ < 9) 

where are phase-space factors for electron and neutrino scattering, the recoil factor R is given 
by Eq.(4), and 07 at is the corresponding elementary cross section for the lepton scattering 
from the moving free nucleon. 

B. Models 

We describe the lepton-nucleon scattering in the impulse approximation (IA), in which 
only one nucleon of the target is involved in the reaction, and the nuclear current is written 
as a sum of single- nucleon currents. Then, the nuclear matrix element in Eq.(7) takes the 
form 

(p,B\J"\A) = J d 3 r exp(rf-r)^ ) (p,r)r At $(r), (10) 
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where T M is the vertex function, t = e B q/W is the recoil-corrected momentum transfer, 
W = a/ (m A + ^>) 2 — Q 2 is the invariant mass, <3> and are relativistic bound-state and 
outgoing wave functions. 

For electron scattering, we use the CC2 electromagnetic vertex function for a free nu- 
cleon [36] 

= F^\Q 2 )Y + ia^^-F^\Q% (11) 

where = i[7 M , 7^/2, Fy 1 ^ and are the Dirac and Pauli nucleon form factors. The 
single- nucleon charged current has V—A structure J M(cc ) = J v + J^. For the free-nucleon 
vertex function T^ cc ^ = T v + Y A we use the CC2 vector current vertex function 

T v = F V (Q 2 )^ + ia^^F M {Q 2 ) (12) 

and the axial current vertex function 

= F A (Q 2 )Yl 5 + F P (Q 2 )q^ 5 . (13) 

The weak vector form factors Fy and Fm can be expressed in terms of the corresponding 
electromagnetic factors for proton F^ and neutron F-^ as follows 

Fi = Fff-Fff. (14) 

For the electromagnetic and weak CC vector vertexes we employ the de Forest prescrip- 
tion [36] (because the bound nucleons are off shell) and Coulomb gauge. For the Dirac 
and Pauli nucleon form factors we use the approximation from Ref. [37] and the dipole 
approximation for the axial Fa and psevdoscalar F P form factors 

Fa{q) - {i+QVMir p{Q) ~ <+Q 2 ' ( 5) 

where F A (0) = 1.267 and m n , Ma are the pion and axial mass, respectively. 

In Ref. [38], a formalism was developed for the A(e, e'N)B reaction that describes channel 
coupling in the FSI of the N + B system. In this work the independent particle shell model 
(IPSM) is assumed for nuclear structure. The model space for 12 C(l,l'N) consists of lsi/2 
and lps/2 nucleon-hole states in the n B and n C nuclei. The l«i/2 state is regarded as a 
discrete state even though its spreading width is actually appreciable. 

In the independent particle shell model the relativistic bound-state functions <3> in Eq.(lO) 
are obtained within the Hartree-Bogolioubov approximation in the a — u model [39]. The 
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upper component of the bound-state wave function $ is used for calculation of the shell 
nucleons spectral function in the PWIA calculations. We use the nucleon bound-state func- 
tions calculated by the TIMORA code [40] with the normalization factors S(a) relative to 
full occupancy of the IPSM orbitals of 12 C: S(lp^/2)=8A%, S , (ls 1 / 2 )=100%, and an average 
factor of about 89%. These estimations of the depletion of hole states follow from the RD- 
WIA analysis of 12 C(e, e'p) for Q 2 < 2 (GeV/c) 2 [41] and consist with a direct measurement 
of the spectral function using 12 C(e, e'p) in parallel kinematics [42], which observed approx- 
imately 0.6 protons in a region with p m > 240 Mev/c and e m > 50 MeV attributable to a 
single-nucleon knockout from correlated cluster. Similar estimates of the depletion of hole 
states are available from the self-consistent Green's function method [43], correlated basis 
function theory [44] and other method also. 

In the RDWIA the ejectile wave function ^ in Eq.(10) is obtained following the direct 
Pauli reduction method [45, 46]. It is well known that the Dirac spinor 



(16) 



can be written in terms of its positive energy component as 

\E+M+S-V +/ 

where S = S(r) and V = V(r) are the scalar and vector potentials for the nucleon with 
energy E. The upper component ^ + can be related to the Schrodinger-like wave function £ 
by the Darwin factor D(r), i.e. 

* + = s/W) t (18) 

E + M + S(r)-V(r) 
K ' E + M V ; 

The two-component wave function £ is the solution of the Schrodinger equation containing 
equivalent central and spin-orbit potentials, which are functions of the scalar and vector 
potentials S and V, and are energy dependent. We use the LEA program [33] for numerical 
calculation of the distorted wave functions with EDAD1 SV relativistic optical potential [47] 
for carbon. 
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In the Plane- Wave Impulse Approximation (PWIA) the final state interaction between 
the outgoing nucleon and the residual nucleus is neglected, and the nonrelativistic PWIA 
exclusive cross section has a factorized form [48] 

1 ^ wr K^ m V(E, p) (20) 

where V(E, p) is the nuclear spectral function. 

According to the JLab data [42, 49], the occupancy of the independent particle shell 
model orbitals of 12 C equals about 89%, on the average. In this work we assume that the 
missing strength (11%) can be attributed to the short-range iVA-correlations in the ground 
state, leading to appearance of high-momentum (HM) and high-energy nucleon distribution 
in the target. In order to estimate this effect in the inclusive cross sections, we consider 
the phenomenological model [34, 35] where the high-momentum (HM) part of the spectral 
function is determined by excited states with one or more nuclei in a continuum. 

We calculate the inclusive cross sections with the FSI effects in the presence of the short- 
range iViV-correlations, using the approach which was proposed in Ref.[21]. 



III. RESULTS AND ANALYSIS 



A. Electron scattering 

The LEA code was successfully tested against 12 C(e,e'p) data. For illustration, Fig.l 
shows measured JLab [49] reduced cross sections for the removal of protons from the Is 
and lp shells of 12 C as functions of missing momentum p m as compared with LEA code 
calculations. It should be noted that negative values of p m correspond to = 7r and positive 
ones to 0=0, where is the angle between the scattering (ki,kf) and reaction {p x ,Pb) 
planes. The data for beam energy Eb eam =2A4:5 GeV and Q 2 =0.6, 1.2, and 1.8 (GeV/c) 2 
were measured in the quasi-perpendicular kinematics with constant (u,q). The detailed 
analysis data [41] for 12 C(e, e'p) with Q 2 <2 (GeV/c) 2 using the RDWIA based upon Dirac- 
Hartree wave functions has shown that the lp normalization extracted from data for Q 2 >0.6 
(GeV/c) 2 is equals approximately 0.87, independent of Q 2 . The total lp and Is strength 
for e m <80 MeV approaches 100% of IPSM, consistent with a continuum contribution for 
30< e m <80 MeV of about 12% of IPSM. 
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FIG. 1: (Color online) Comparison of the RDWIA and the RFGM calculations for electron, neu- 
trino and antineutrino reduced cross sections for the removal of nucleons from lp and Is shells 
of 12 C as functions of missing momentum. JLab data [49] for beam energy Eb eam =2A55 GeV, 
proton kinetic energy T p =350 MeV, and Q 2 =0.64 (GeV/c) 2 . The RDWIA calculations are shown 
for electron scattering (dashed line) and neutrino (solid line) and antineutrino (dashed-dotted line) 
scattering; and the RFGM results are shown for the reduced cross sections (dotted line) for the 
JLab kinematics. 

The electron and neutrino scattering off the nuclei are closely interrelated and one can 
treat both processes within the same formalism. In the nonrelativistic PWIA, a re d is a nu- 
clear spectral function and should be similar for electron and (anti)neutrino scattering except 
small distinctions which can be attributed to the Coulomb distortion upon the electron wave 
function. The small difference between neutrino and antineutrino is due to difference in the 
FSI of the proton and neutron with the residual nucleus. This effect is neglected at the 
energy beam higher than 1 GeV. There is an overall good agreement between calculated in 
the RDWIA electron and (anti)neutrino cross sections and data. Apparently the RFGM 
predictions (with the Fermi momentum pp=221 MeV/c and binding energy 6^=25 MeV) 
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overestimate the values of cross sections and completely off the exclusive data. This is 
due to the uniform momentum distribution of the Fermi gas model and neglecting by the 
FSI effects. Therefore, the RFGM can not predict well enough the momentum distribution 
of outgoing protons in simulation of the CC QE two-track events at momentum transfer 
\q\ < \p m \, i.e. at low Q 2 . 

A complex relativistic optical potential with a nonzero imaginary part generally produces 
an absorption of flux. However, for the inclusive reaction, the total flux must conserve. 
Currently there is no fully consistent solution to this problem, and different approaches are 
used. In Refs.[15, 50] it was shown that the inclusive CC cross sections calculated with 
only the real part of the optical potential are almost identical those of the Green's function 
approach, in which the FSI effect in the inclusive reactions is treated by means of a complex 
optical potential and the total flux is conserved. In this work, in order to calculate the 
inclusive and total cross sections, we use the approach, in which only the real part of the 
optical potential EDAD1 is included. 

To test our approach, we calculated the inclusive 12 C(e, e') cross sections and compared 
them with data from SLAC [51-53], from Saclay [54], and from JLab [55]. Figures 2 
and 3 show measured inclusive cross sections as functions of energy transfer as compared 
to the RDWIA, PWIA, and RFGM calculations. These data cover the range of the tree- 
momentum transfer (around the peak) form \q\ «310 MeV/c (beam energy E e =1500 MeV 
and scattering angle 0=11.95°) up to \q\ ^530 MeV/c (£ e =2020 MeV, 0=15°). We note 
that, relative to the PWIA results, the generic effect of the FSI with the real part of the 
optical potential is to reduce the cross section value around the peak and to shift the peak 
toward the lower value energy transfer. The inclusion of the high-momentum component 
increases the inclusive cross section in the high-energy transfer region and improves the 
agreement with data. For the RDWIA results, the difference between the calculated and 
measured cross sections at the maximum is less than ±12%. For the RFGM results these 
difference decreases with \q\ from about 20% at \q\ ?a310 MeV/c down to «13% at \q\ ?a510 
MeV/c. These results demonstrate a strong nuclear-model dependence of the inclusive cross 
sections at low momentum transfer. 
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FIG. 2: (Color online) Inclusive cross section versus energy transfer u for electron scattering on 12 C. 
The data are from Ref.[51] (filled triangles), Ref.[52] (open squares), and Ref.[54] (filled squares). 
In Ref.[51] data are for the electron beam energy £^=1500 MeV, and scattering angles # e =11.95°, 
13.54°; in Ref.[52] data are for £ e =537 MeV and 6» e =37.1°; in Ref.[54] data are for £ e =361 MeV 
and 9 e =60°. 
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FIG. 3: (Color online) Same as Fig. 2, but the data are from Ref.[54] (filled squares) for £"6=620 
MeV, (9 e =36° and £ e =680 MeV, (9 e =36°; Ref.[53] (filled circles) for £ e =500 MeV, 6> e =60°; Ref.[55] 
(stars) for £ e =2020 MeV, e =15°. 

B. Neutrino scattering 

The charged-current QE events distributions as functions of Q 2 were measured in K2K [12, 
56] and MiniBooNE [1] experiments. The shape of the Q 2 distribution, which is weakly 
dependents on the flux uncertainties, was analyzed. High statistic data show a disagreement 
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with the RFGM predictions. The data samples exhibit significant deficit in the region of 
low Q 2 <0.2 (GeV/c) 2 (so call low-Q 2 problem). In Ref.fl] it was shown that the data/MC 
disagreement is not due to mis-modeling of the incoming neutrino energy spectrum, but to 
inaccuracy in the simulation of CC QE interactions. To tune the Fermi gas model to the low 
Q 2 , an additional parameter k was introduced which reduced the phase space volume of the 
nucleon Fermi gas at low-momentum transfer. This parameter controls the Q 2 distribution 
in the low Q 2 region only. 

In the region of high Q 2 the data excess is observed, and the value of the axial vector 
mass Ma obtained from a fit to the measured data, are higher than the results of previous 
experiments. The formal averaging of Ma values from several experiments, which are very 
wide spread from 0.7 to 1.3 GeV, was done in Ref.[57]: M J 4=1.026±0.021. This result is 
also known as the axial mass world average value. K2K obtained the value of 1.2±0.12 
from the SciFi detector [12] using the water-aluminum mixture as a target, and also the 
preliminary result 1.14±0.11 from the SciBar detector [56] using a scintillator target. The 
MiniBooNE experiment (scintillator target) found that the data were better described with 
an adjustment of two parameters M^=1.23±0.20 GeV and /t=1.019±0.011 [1]. 

Recently the NOMAD experiment [58] extracted the value of M A =1.05±0.02±0.06 GeV 
using a carbon target, which is in agreement with the world average value. This result was 
obtained from the analysis of a measured u 12 C total CC QE cross section for neutrino energy 
above ^4 GeV, where the cross section 'plateaus' is reasonably well-known. It should be 
noted that both approaches, i.e. analysis of the shape of the Q 2 distribution and the direct 
measurement of the total cross section assume, that the vector form factors are known well 
from the electron scattering experiments. Actually, at Q 2 > 3 (GeV/c) 2 the values of the 
neutron form factors are much less known that those of the proton [59], and the relative 
contribution from this region to the total cross section increases with neutrino energy. 

To study nuclear effects on the (^-distribution, we calculated with M A = 1.032 GeV the 
inclusive cross sections da/dQ 2 for neutrino energies £^=0.5, 0.7, 1.2 and 2.5 GeV and 
compared them with those for neutrino scattering on a free nucleon. The results for neutrino 
and antineutrino scattering on carbon are presented in Figs. 4 and 5, respectively, which 
show da/dQ 2 as functions of Q 2 . Here, the results obtained in the RDWIA, are compared 
with cross sections calculated in the PWIA and RFGM. The cross sections for the exclusive 
reaction are shown as well. In the region Q 2 < 0.2 (GeV/c) 2 the Fermi gas model results 
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FIG. 4: (Color online) Inclusive cross section vs the four-momentum transfer Q 2 for neutrino 
scattering off 12 C and for the four values of incoming neutrino energy: £^=0.5,0.7,1.2, and 2.5 
GeV. The solid line is the RDWIA calculation, whereas the dashed and dash-dotted lines are, 
respectively, the RFGM and PWIA calculations. The dotted lines are the cross sections for the 
exclusive reaction. 
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FIG. 5: (Color online) Same as Fig. 4, but for antineutrino scattering. 

for neutrino (antineutrino) are higher than those obtained within the RDWIA. At Q 2 =0.1 
(GeV/c) 2 this discrepancy equals 12% (28%) for £^=0.5 GeV and decreases to 7% (12%) for 
£j,=2.5 GeV. The contribution of (v, nN) channels to the inclusive cross sections is about 
60%. 

Nuclear effects on the shape of the four-momentum transfer Q 2 distribution, i.e. ra- 
tio R(e u ,Q 2 ) = (da/dQ 2 ) nuc /(da/dQ 2 )f ree , where (da/dQ 2 ) nuc is the cross section scaled 
with number of neutron/proton in the target and (da / dQ 2 ) f ree is the cross section for 
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FIG. 6: (Color online) Ratio R[e v ,Q ) vs the four-momentum transfer Q 2 for neutrino scattering 
off 12 C and for the four values of incoming neutrino energy: £^=0.5,0.7,1.2, and 2.5 GeV. As shown 
in the key, the cross sections were calculated with the RDWIA, PWIA, and RFGM. 

(anti) neutrino scattering off free nucleon, are presented in Fig. 6 for neutrino and in Fig. 7 
for antineutrino as functions of Q 2 . Here, the results obtained in the RDWIA for energies 
e u =0.5, 0.7, 1.2 and 2.5 GeV are compared with those calculated in the PWIA and Fermi 
gas model. The nuclear effects are seen at low Q 2 ; the tail of the momentum distribution 
at high Q 2 , an overall suppression, and a slight change in the slope in the middle region at 
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FIG. 7: (Color online) Same as Fig. 7, but for antineutrino scattering. 

e v >1 GeV can also be seen. The range of Q 2 where R m const, i.e. nuclear effects are 
negligible and, therefore cannot modify the value of Ma, increases with incoming neutrino 
energy. At energy higher than 1 GeV the range 0.3 < Q 2 < 1 (GeV/c) 2 can be used for M4 
extraction from Q 2 shape-only fit. 

We calculated da/dQ 2 cross sections at energy 700 MeV in the RDWIA and Fermi gas 
model with Ma= 1.032 and 1.32 GeV. The results are shown in Fig. 8 as functions of Q 2 . 
Apparently at low Q 2 <0.1 (GeV/c) 2 the cross sections depend weakly on the value of the 
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FIG. 8: (Color online) Inclusive cross section vs the four-momentum transfer Q 2 for neutrino 
(upper panel) and antineutrino (lower panel) scattering off 12 C with energy e v =0.7 GeV and for 
the two values of axial mass Ma= 1.032 and 1.32 GeV. As shown in the key, cross sections were 
calculated within the RDWIA and Fermi gas model. 

axial mass and Q 2 distributions are controlled by nuclear effects. 

The inclusive neutrino and antineutrino cross sections for energies e u =0.5, 0.7, 1.2, and 
2.5 GeV are presented in Figs. 9 and 10, which show da/de^ as a function of muon energy. 
Here, the results obtained in the RDWIA with Ma= 1.032 GeV are compared with the 
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FIG. 9: (Color online) Inclusive cross section vs the muon energy for neutrino scattering on 12 C 
and for the four values of incoming neutrino energy: £^=0.5,0.7,1.2, and 2.5 GeV. As shown in 
the key, the cross sections were calculated with the RDWIA, PWIA, RFGM, and RDWIA for the 
exclusive reaction. 

inclusive cross sections calculated in the PWIA, RFGM, and RDWIA for the exclusive 
reaction. The cross section values obtained in the RFGM are higher than those obtained 
within the RDWIA. For neutrino (antineutrino) cross sections in the region close to the 
maximum this discrepancy is about 25%(49%) for £^=0.5 GeV and 23%(29%) for e v =2.5 
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FIG. 10: (Color online) Same as Fig. 9, but for antineutrino. 

GeV. The contribution of (u, fiN) channels to the inclusive cross sections is about 60%. 

The neutrino and antineutrino total cross sections calculated with M^= 1.032 GeV up 
to neutrino energy 2.8 GeV, are shown in Fig. 11 together with data of Refs. [60-63] . Also 
shown are the results obtained in the RFGM, PWIA as well as the contribution of the 
exclusive channels to the total cross sections. The cross sections are scaled with the number 
of neutron/proton in the target. The ratio between the neutrino cross sections calculated in 
the RFGM and RDWIA decreases with neutrino energy from about 1.15 for e u =0.5 GeV to 
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FIG. 11: (Color online) Total cross section for CC QE scattering of muon neutrino (upper panel) 
and antineutrino (lower panel) on 12 C as a function of incoming (anti) neutrino energy. The solid 
line is the RDWIA result while the dashed-dotted and dotted lines are, respectively, the RFGM 
and PWIA results. The dashed line is the RDWIA result for exclusive reaction. Data points for 
different targets are from Refs. [60-63] 
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FIG. 12: (Color online) Ratio of the total cross sections per neutron/proton R=0/C for CC QE 
scattering of muon neutrino (upper panel) and antineutrino (lower panel) scattering on 16 and 
12 C vs incoming (anti)neutrino energy. The solid line is the RDWIA result, while the dashed 
and dashed-dotted lines are, respectively, the RFGM and RDWIA without contributions of the 
short-range correlations. 

1.02 for e„=2.6 GeV. For the antineutrino cross sections this ratio is about 1.5 for £^=0.5 
GeV, and 1.05 for £^=2.6 GeV. The contribution of the exclusive channels is about 60%. 
The results presented in Fig. 11 show significant nuclear-model dependence for energy less 
than 1 GeV. 

The RDWIA prediction for the CC QE flux-averaged total cross section is compared with 
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the experimental result from the LSND Collaboration at the Los Alamos for 12 C(z/ M , /x ) re- 
action [64]. The mean energy of the neutrino flux above threshold is 156 MeV. The calculated 
value of 10.14xl(r 40 cm 2 well agree with measured value of (10.46±0.3±1.8)xl0~ 40 cm 2 . 

To compare the CC QE total cross sections for (anti) neutrino scattering on the oxygen [21] 
and carbon targets, we calculated the ratio R(e„) = {<Jtot)nud/{&tot)nuci-> where the cross 
sections {cr l tot ) nuc i are scaled with the number of neutron/proton in the targets. The results 
obtained in the RFGM and RDWIA are shown in Fig. 12. The Fermi gas model predicts 
almost identical values of (al ot ) nuc i for 16 and 12 C. In the RDWIA approach the cross 
section calculated for oxygen is lower than that for carbon. For the neutrino (antineutrino) 
scattering this ratio is 0.90(0.88) at £^=0.7 GeV and 0.92(0.91) at £„=2.6 GeV. 

To study the iViV-correlation effects, we calculated the ratio R(e u ) without the NN- 
correlation contribution, i.e. with S a —1 for all bound nucleon states in the oxygen and car- 
bon targets. The difference between results obtained with and without the high-momentum 
component contribution is about 5% for e v >1 GeV. In Ref.[21] was shown that the NN- 
correlation effect reduces the total cross section in proportion to the missing strength in the 
nuclear ground state, which is about 25% for 16 and 11% for 12 C. 

Therefore in the long-base line neutrino oscillation experiments a part of the near detector 
must include some of the same target material, as the far detector to reduce the systematic 
uncertainty due to nuclear effects on the CC QE total cross section. 

IV. CONCLUSIONS 

In this paper, we study electron and CC quasi-elastic (anti)neutrino scattering on a carbon 
target in different approximations (PWIA, RDWIA, RFGM) placing particular emphasis 
on the nuclear-model dependence of the results. In RDWIA, the LEA program, adapted 
to neutrino interactions, was used to calculate the differential and reduced exclusive cross 
sections. We found that the reduced cross sections for (anti) neutrino scattering are similar 
to those of electron scattering and the latter are in a good agreement with electron data. 
In calculating the inclusive and total cross sections, the imaginary part of a relativistic 
optical potential was neglected and the SRC effect in the target ground state was taken into 
account. This approach was tested against electron-carbon inclusive scattering data. This 
test revealed an overall agreement with the data, with the differences between calculated 
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and measured cross sections in the peak region less than 12%. 

We calculated da/dQ 2 cross sections for different neutrino energies and estimated the 
range of Q 2 where nuclear effects on the shape of Q 2 distribution are negligible. Also was 
shown that at low Q 2 <0.1 (GeV/c) 2 the coss sections depend weakly of the values of the 
axial mass. 

The CC QE total cross sections predicted by the RFGM are higher than the corresponding 
values obtained in the RDWIA and this difference decreases with neutrino energy. The flux- 
averaged total cross section was calculated within the RDWIA approach and compared with 
the experimental result from the LSND Collaboration. The calculated cross section is in 
good agreement with data. We compared the CC QE total cross sections (scaled with the 
number of neutron/proton in the target) for (anti) neutrino scattering on the oxygen and 
carbon targets and found that the cross sections calculated within the RDWIA for oxygen, 
are lower than those calculated for carbon, and the SRC effects increase this difference. 

We conclude that the data favor the RDWIA results. This indicates that the use of 
RDWIA in Monte Carlo simulations of the neutrino detector response would allow one to 
reduce the systematic uncertainty in neutrino oscillation parameters. 
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